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ABSTRACT 

ERTS-1 multispectral scanner (MSS) digital data was used to 
delineate coastal land, vegetative, and water features in two 
portions of the Texas Coastal Zone. ERTS-1 MSS data (Scene ID's 
1037-16244 and 1037-16251) acquired on August 29, 1972, were 
analyzed on NASA Johnson Space Center systems through the use of 
two clustering algorithms. Seventeen to 30 spectrally homogeneous 
classes were so defined. Many classes were identified as being 
pure features such as water masses, salt marsh, beaches, pine, 
hardwoods, and exposed soil or construction materials. Most 
classes were identified to be mixtures of the pure class types. 
Using an objective technique for measuring the percentage of 
wetland along salt marsh boundaries, an analysis was made of the 
accuracy of areal measurement of salt marshes. Accuracies ranged 
from 89 to 99 percent. .NASA aircraft photography was used as the 
basis for determining the true areal size of salt marshes in the 
study sites. 


1. INTRODUCTION 

The purpose of this paper is to report on the application of two classifi- 
cation algorithms to ERTS-1 multispectral scanner (MSS) system corrected data 
acquired over the Galveston Bay (Texas) area on August 29, 1972 (Scene ID's 
1037-16244 and 1037-16251). Specifically, two algorithms (ISOCLS and NSCLAS) 
available at the NASA Johnson Space Center (JSC) were used to produce classi- 
fication symbol maps of portions of the Texas coastal zone. The resulting 
classes then were related to an information hierarchy selected for the MSC 
ERTS-1 coastal estuarine investigation. 
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Nonwater 
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Wetland 

Other 


Level 3 
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Nonproductive water 

Vegetated 

Nonvegetated 
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Photography obtained by the NASA WB-57 aircraft on August 30, 1972, was 
used to aid in the identification of the classes defined by the algorithms 
and to provide a basis for the measurement of the true area of selected 
features such as water bodies and marsh. 
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2. METHOD OF ANALYSIS 

Two study sites were edited from ERTS-1 MSS system corrected data contained 
in frames 1037-16244 and 1037-16251. These sites covered (1) Trinity River 
Delta area including salt, brackish, and fresh water marshes, and (2) Galveston 
Island, West Bay, Pelican Island, and offshore Gulf waters. 

One computer compatible tape (CCT) containing the Galveston data was sent 
to Purdue University. This data set could then be analyzed through the use of 
the remote terminal at JSC. Both sets could be analyzed by JSC systems. 

Since the two study sites exhibited natural features of complex configura- 
tion and unknown spatial extent, it was decided to use classification algorithms 
to produce classification symbol maps of the sites. Each of the two algorithms 
(ISOCLS and NSCLAS) attempts to partition the data into spectrally homogeneous 
groups. They then generate classification maps to display the geographical 
regions in which these homogeneous groups of data occur. A detailed descrip- 
tion of ISOCLS can be found in references 1 through 5. A description of NSCLAS 
is given in reference 6. 

Photography (film 2443, color IR) of the study sites was obtained by the 
NASA WB- 57 aircraft from 60,000 feet (MSL) on August 30, 1972, one day after 
the subject ERTS-1 data. The photography was rectified and enlarged to a 
scale of 1:41,200. Then, wetland areas were delineated by coastal analysis 
team members. The area of each delineation was determined through the use of 
a photographic data quantizer (planimeter) . Field trips into the sites were 
made to aid in the identification of classes defined by ISOCLS in that site. 
Photo interpreters identified classes in both sites through the use of the 
aircraft photography. 

In the case of the areal measurement of the wetland areas in the Galveston 
site, an objective technique described in reference 7 in this symposium was 
used to count picture elements and partial picture elements on the boundary of 
the marsh. 


3 . RESULTS 

Space does not permit the inclusion of all of the color coded classifica- 
tion maps produced during the present analysis effort. Examples of level two 
maps in each site are shown in figures 1 and 2 (originals are in color). In 
figure 2, four classes are shown: water, forested areas, vegetation including 

marsh, and bare soil/beach/asphalt/concrete. In figure 3, three classes are 
shown: water, wetland (marsh), and other. 

4. WETLAND ANALYSIS 

Using the objective technique described in reference 7 (this symposium), 
the areas of three wetlands were determined. The results are given in table I. 

TABLE I.- RESULTS OF DETERMINING AREA OF WETLANDS FROM 
ERTS-1 DATA THROUGH AN OBJECTIVE TECHNIQUE 


ESTIMATED AREA PERCENT ACCURACY 


WETLAND 

NUMBER 

(ACRES) 

ISOCLS NSCLAS 

TRUE AREA 
(ACRES) 

ISOCLS 

NSCLAS 

1 

864.13 

871.56 

866.0 

99.6 

99.4 

2 

197.97 

193.77 

191.3 

96.5 

98.7 

3 

150.86 

165.94 

149.6 

99.2 

89.1 

The 

accuracies in wetland (marsh) 

area measurement 

are quite high. 

Both 


algorithms accurately define wetland features and boundary classes. It is not 
possible at this time to determine the more superior algorithm; however, 

ISOCLS on the UNIVAC 1108 was able to handle a larger number of data points than 
NSCLAS on the IBM 360/67 at Purdue University. 

5. MIXTURE PROBLEM 

A graph depicting the position of cluster means in the first and fourth 
channels of the ERTS-1 MSS is shown for the Trinity site in figure 4 for a 
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30-cluster case, which will be the subject for the remainder of this discussion. 
This figure is a good representation of the data and of the results of the 
cluster analysis. With few exceptions the second and third channels of the 
MSS add little new information. 

The clusters obtained here are, for the most part, mixtures of some of the 
major classes that exist on the surface. The remaining clusters can be uniquely 
related to these major classes. The following table shows those clusters that 
can be so related: 

Cluster Identification Class 

B,R,D,N,S Vegetation excluding trees 

C Pine Trees 

4 Hardwoods 

F Beach/asphalt/concrete 

U,K Bare soil 

1 ,G,H,L,S Water 

The clusters have been identified by assigning to them consecutive arbitrary 
symbols 1, 2, .... 9, A, B, ... 

Up to five clusters can be identified as belonging to water bodies. No 
quantitative relationship can be assigned to these clusters at this time 
because there was no concurrent water survey on the date of the ERTS pass. 

The following is a general description of the type of water clusters and where 
they occur : 

Cluster Identification 

1 

G 
H 


L 

S 


Class or Description 

Bay water (probably slightly 
turbid) 

Turbid (fresh or saline) 
water 

Highly turbid water, very 
shallow water, or under- 
water spoil banks 

Fresh water that occurs in 
lakes 

Turbid fresh water or turbid 
water with a small amount 
of vegetation (less than 
10 percent) . 


The clustering analysis differentiates forested areas from other types of 
vegetation. Within the forested areas, pine trees and hardwoods can also be 
differentiated and assigned to separate clusters. Vegetation types other than 
trees cannot be separated. Specifically, the same combination of clusters 
occur over vegetated areas which are known to contain different classes such 
as Phragmites communes . Spartina spartinae , and pasture land. 

One conspicuous fact that this investigation has brought forth is that a 
good part of the data can be related to a mixture of two or more classes found 
on the surface. In the case of wetlands area, the proportion of classes con- 
tributing to a particular mixture varies in a continuous fashion over the 
possible range. Under these circumstances there are no natural groupings of 
the data; so the clustering program tends to partition the data at arbitrary 
places . 

It may be noted that data assigned to clusters which represent a mixture 
of classes, can be obtained by linearly combining data assigned to clusters 
identified with one of the major classes. This fact may be easily seen by 
inspecting figure 4 where, for example, data points in cluster 8 can be 
obtained by linearly combining data points in clusters 5 and C. In this manner, 
a range of ratios of a particular mixture can be extracted from the data. The 
following table represents this range for the water/vegetation mixture. These 
numbers represent an estimate only and are not the result of a statistical 
s tudy . 
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Cluster Identification 


Ratio of Water 


6 

Q 

3 

E 

2 

7 

5 

N 

R 


. 6-. 95 
.5-. 8 
.3-. 7 
. 2-. 6 
. 1- . 5 
0-.35 
0-.20 
0-.10 
0-.05 


Other mixtures could consist of almost any possible combination of the 
major classes that have been identified. The following table presents a list 
of possible mixutres that have been verified: 


Cluster 

Identification 


Mixture 


4 

C 

0 

8 


6,Q,3,2,7 

1,9 

P, A, T, M, J 
E,U,K 


Hardwoods; pines and vegetation (excl. trees) 

Pines; hardwoods and water 
Dense pines; trees and water 

Pines and vegetation (excl. trees); hardwoods vegetation 
(excl. water) and water 
Marsh; vegetation and water 

Water and bare soil/beach/asphalt in increasing proportion of 
bare soil/beach/asphalt 

Vegetation and bare soil/beach/asphalt/concrete in increasing 
proportion of bare soil/beach/asphalt/concrete 
Vegetation, water, and bare soil/beach/ asphalt/concrete in 
increasing proportion of bare soil/beach/asphalt/concrete 


6. CONCLUSION 

Two clustering algorithms were used to produce classification maps of the 
Trinity River delta area (Scene ID 1037-16244) and the Galveston area (scene 
ID 1037-16251) from data obtained on August 29, 1972. Seventeen to thirty 
classes were defined. Most of these classes represented mixtures of land, 
water, and vegetation including forest. Through the use of a cluster diagram, 
ranges of percent of water and wetland were assigned to some classes. Improve- 
ments would result if a statistical approach were used for the mixture problem. 


7. RECOMMENDATIONS 

Further analysis is needed to determine whether the position of cluster 
centers in the cluster diagram can be used to identify the cluster as repre- 
senting a certain class of coastal features. 
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